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Abstract 

Inspired by Holstein's work on small polaron hopping, the evolution equations of localized states 
and extended states in presence of atomic vibrations are derived for an amorphous semiconductor. 
The transition probabilities are obtained for four types of transitions: from one localized state to 
another localized state, from a localized state to an extended state, from an extended state to a 
localized state, and from one extended state to another extended state. At a temperature not too 
low, any process involving localized state is activated. The computed mobility of the transitions 
between localized states agrees with the observed 'hopping mobility'. We suggest that the observed 
'drift mobility' originates from the transitions from localized states to extended states. Analysis 
of the transition probability from an extended state to a localized state suggests that there exists 
a short-lifetime belt of extended states inside conduction band or valence band. It agrees with 
the fact that photoluminescence lifetime decreases with frequency in a-Si/Si02 quantum well while 
photoluminescence lifetime is not sensitive to frequency in c-Si/Si02 structure. 
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I. INTRODUCTION 



From the energy spectrum of an amorphous semiconductor-, one knows that there are 
four types of carrier transitions that may contribute to the electric conduction: type (1) 
transition between two locahzed states; type (2) transition from a locahzed state to an 
extended state; type (3) transition from an extended state to a locahzed state and type (4) 
transition between two extended states. Type (4) transition is common to both crystalline 
and non-crystalline materials. For amorphous semiconductors, type (1) transition is the main 
conduction mechanism in a broad temperature range. It has been investigated in various 
ways^'^. Although the thermal equilibrium population of the extended states is lower than 
the population of the localized states, the contribution from the carriers in the extended 
states to the electric conduction is still observable for a moderately high temperature. An 
investigation of the transitions of type (2), type (3) and type (4) is necessary for a complete 
description of carrier dynamics. 

In other electronic hopping processes, thermal vibrations of atoms also play an important 
role. Electron transfer (in polar solvent and inside large molecules}^ and polaron diffusion in 
a molecular crystal^*^ are two examples. The transition probability Wll between two sites 
in a thermally activated process is given by the Marcus formula 



where ull has the dimension of frequency for a specific hopping process. is the tem- 
perature dependent activation energy. Xll is the reorganization energy, AG^ll is the energy 
difference between the final state and the initial state^. Eq.([T]) is valid for both electron 
transfer- and small polaron hopping-. The mathematical form of Holstein's work for one 
dimensional molecular crystal is quite flexible and can be used to treat three dimensional 
materials with slight modifications^'^iiiSi'iiii^. The effect of static disorder may be taken into 
account by replacing a fixed transfer integral with a distribution. The static disorder reduces 
the strength of the electron-phonon (e-ph) coupling needed to stabilize global small-polaron 
fomation^^. Amorphous semiconductors offer a different regime, in which the static disorder 
is so strong that the states in valence and conduction tails are localized. These localized 
states interact with the atomic vibrations. 

We will extend Holstein's work-i^ to amorphous semiconductors. In Sec. II, we first intro- 
duce some notation about localized states, extended states and electron-phonon interaction. 




(1) 
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Then the equations of time evolution for locahzed states and extended states in presence 
of atomic vibrations are derived from a time-dependent Schrodinger equation. The con- 
nections with electron transfer, with Kramers' problem of escape across a barrier and with 
small polar ons are pointed out. 

In Sec. Ill we study the hopping processes among localized states (LL) in amor- 
phous solids. Eq.([T]) is re-established in high temperature regime (A;_bT > 2.5JtIj, ZJ = 
J dujD{uj)uj/ J dujD{uj) is the average phonon frequency, D{uj) is the spectral distribution 
of phonons, the factor 2.5 comes from the requirement that the error of approximation 
csch^l^ « is l^ss than 0.003, cf. the paragraph below Eq.( l84l) ). The reorganization 
energy \ll is expressed by the phonon spectrum, eigenvectors of the normal modes and 
the electron-phonon interaction parameters. The computed temperature dependence of the 
mobility of LL transition in a-Si agrees with that of the observed 'hopping' mobility for two 
regimes T>250K and T<250K (cf. Fig {Hand Figl2]). 

The transition from a localized state to an extended state (LE) induced by the transfer 
integral is reported in Sec. IV, it has been suggested as the main conduction mechanism in 
amorphous silicon and is called 'phonon induced delocalization'^^. Without the dressing 
effects of the vibrations of atoms, in certain sense, type (1) transition is similar to the 
transition between two bound states in a molecule and to the electron transfer between two 
ions, type (2) transition is analogous to the ionization process of an atom and to the escape 
across a potential barrier. 

The transition from an extended state to a localized state (EL) induced by the electron- 
phonon interaction, is presented in Sec.V. In conduction band, the energy of an extended 
state is higher than that of a localized state in the lower tail. For an extended state with 
high enough energy, EL transition is in Marcus inverted regime. One expects that there 
exists a short-lifetime belt of the extended states inside conduction band or valence band 
(cf. FiglH]). The states inside this belt favor non-radiative transitions by emitting several 
phonons. The conjectured short-lifetime belt in conduction band agrees the fact that the 
photoluminescence lifetime of a-Si/Si02 quantum well decreases with frequency while that 
of c-Si/Si02 is not sensitive to frequency. Type (3) transition is similar to the free electron 
capture process in an atom and to the capture a particle by a well in a viscous liquid. 
One common feature of the type (1), (2) and (3) transitions is that at higher temperature 
ksT > 2.5hLJ, the probability of any transition involving localized state takes the form of 
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Eq.([T]). We suggest that the observed drift mobihty may originate from type (2) transitions. 
Sec. VI is devoted to the transition between two extended states caused by electron-phonon 
interaction. Higher order processes and conduction mechanisms are briefly discussed. Finally 
we summarize this work and mention some un-touched problems. 

II. EVOLUTION OF THE STATES DRIVEN BY THE VIBRATIONS 
A. Localized and extended states 

Consider an amorphous sample with M atoms. Denote the static positions of the atoms 
as {T^n; n = 1;2, ■■■ jA/"}, the displacements of the atoms due to thermal vibrations as 
{u", n = 1, 2, ■ ■ ■ , A/"}. To make the formulae compact, we rename the SA/" vibrational de- 
grees of freedom [u]., Uy,ul,--- , u^, u^, } as {xk, k = 1,2, ■ ■ ■ , 3Af} and rename the 3Af 
static position coordinates {TZix, T^iy, T^iz, ■ ■ ■ , T^aTx, T^Afy, T^aTz} as {Xj, j = 1,2, ■ ■ ■ , 3J\f}. 
Consider the single electron Hamiltonian: 

^« = -|^v^ + E^^(r-7^n) (2) 

n=l 

where U{r — TZ^) is the potential energy felt by an electron at r due to an atom at ^^y^. 
In Eq.([2]) only the static disorder of the amorphous sample is taken into account, ha is a 
mean field approximation, the existence of U is known from Hartree-Fock method or density 
functional theory, ha has two kinds of eigenstates: localized states {0^^} 

ha<PX (r, {7^„}) = El ({7^„})0° ^ (r, {7^„}) (3) 

and extended states {C,Bi} 

haU (r, {7^„}) = Eb-, ({7^n})U (r, {7^„}) (4) 

Eigenstates belonging to different eigenvalues are orthogonal to each other. In conduction 
band, the energies of the extended states are above those of the localized states. We use 
Ai,A2--- = 1, 2, ■ ■ ■ , Nl to label the localized states (A'^^: is the total number of the localized 
states), i?2 ■ ■ ■ = 1, 2, ■ ■ ■ , A''^; to label the extended states {Ne is the total number of 
extended states). 

When an electron in an extended state is scattered by phonons, the resulting state is still 
an extended state. The extended states are not radically modified by thermal vibrations. 
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The situation for the locahzed states is different. The wave function of a locahzed state 
is non-zero only in some finite spatial region Da^'- only the vibrations of the atoms inside 
Da^ couple to 0^^: 

I dvct^%[u-.VJJ{v-n^)]ctPA, (5) 

The 'concentrated' wave function makes ([5]) comparable to or even larger than the transfer 
integral between two localized states and the transfer integral between a localized state and 
an extended state. The change in wave function of a localized state must be taken into 
account in the zeroth order: 

[-|^^'+ E u{r-n^,vi-)]<PAA^-nA,,{x^i^}) = EAA{^X'^)<pM^^ (6) 



he 

n 



where TZai is a fixed point (could be arbitrarily chosen) inside -D^i, Pa^ is the index of t 
vibrational degrees of freedom inside Da^- Eq.([n]) is a generalization of Eq.(I-3) in ref|5|. 
The effect of the atoms in the neighboring regions of Da^ is neglected. The displacements 
of the atoms change the overlap integrals between localized states 

Ta,a, = j dhct>X{r - 7^^„ {x^jj)0^,(r - 7^^„ = 5a,a, + Sa,a, (7) 

For two localized states A2 and Ai, SA2A1 is a small quantity if Da^ and Da2 do not overlap. 
Using the linear approximation for e-ph interation^, 



(8) 

Eq.([8]) is the first order correction of electronic energy due to the e-ph interaction. The 
most localized states can be described by ,^~'^/^e~^''^, where ^ is localization length of a 
localized state, r is the distance between electron and some representative point inside -Dai- 
dpA^ ~ 4^*Qg2 ; Co is the permittivity of vacuum, Z* is the effective nuclear charge of an atom 
felt by a conduction electron. For electrons in a metal, researchers usually focus on how an 
electron in an extended state is scattered into another extended state by the e-ph interaction 
rather than the correction to energy. Unlike , includes the effect of vibrations of the 
atoms and is no longer orthogonal to any extended states of ([2]), 

j d'rCn^ir, {7^„})0A.(r - 7^A„ = Y^.a, (9) 
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B. Single-electron approximation 



For definiteness, we consider the electrons in the conduction band of an amorphous semi- 
conductor. For carriers in the mid-gap states and the holes in valence band, we need only 
slightly modify the notation. In intrinsic and lightly doped n-type semiconductors, the num- 
ber of the electrons is much smaller than the number of the localized states. The correlation 
between electrons in a hopping process and the screen effect caused by these electrons can 
be neglected. Essentially we have a single particle problem: one electron moves in many 
empty localized states and extended states in the conduction band. 

Consider one electron moving in an amorphous solid with M atoms, the total Hamiltonian 
of the system is 

^=^^' + i;^(r'Wn) + E^V^+ E nW„,W„0 (10) 

n=l n " n,n'(n7^n') 

where n denote the site of network, Wn is the instantaneous position vector of the nucleus 
at site n, U (r, Wn) is the interaction between the electron and the nucleus at the n**^ site. 
^(Wn, Wn') is the total effective interaction of the nucleus at site n' and the nucleus at site 
n, which including both the Coulomb repulsion between them and the induced attraction 
by the electrons. The total wave function ?/'(r, {Wn};t) of the system is a function of all 
degrees of freedom, its time evolution is determined by 

z;i^^(r, {Wn}; t) = Hi^ir, {Wn}; t). (11) 

For temperature well below melting point, the system executes small harmonic oscillations. 
The motion of nuclei Wn = TZ^ + u" can be viewed as vibrations of the atoms {u"} around 
their equilibrium positions {TZn}- ^(r, {Wn}; t) is changed into ip{r, {u"}; t), a function of 
displacements of the atoms. H is simplified as 



Hi = he + K (12) 



where 

he 



n=l j ^ jk 

hf. is the single electron Hamiltonian including the vibrations, ha in Eq.([2]) and the Hamil- 
tonian in Eq.([6]) are two different approximations of h^- hy is the vibrational Hamiltonian, 



(kjk) is the matrix of force constants. The evolution of the total wave function -^(r, {u"}; t) 
of the system of "one electron + many nuclei" is given by Schrodinger equation 

«^^^^^ = //,*(r. {«"};*) (14) 



C. Evolution equations 

The Hilbert space of is spanned by the localized states and the extended states. The 
total wave function ?/'(r, {u"}; t) of the system of "one electron + many nuclei" can be 
expanded as 



V^(r,xi, ■ ■ ■ ,X3j^;t) = ^aAi(xi, ■ ■ ■ ,X3Ar;t)0Ai + ■ ■ ■ ,X3j^;t)^Bi (15) 

where a^^ is the probability amplitude at moment t that the electron is in localized state 
Ai while the displacements of the nuclei are {xj, j = 1,2, ■ ■ - SAf}, 6^1 is the amplitude at 
moment t that the electron is in extended state Bi while the displacements of the nuclei are 
{xj, j = 1,2, ■ ■ ■ 3Af}. The first sum runs over all the localized states, the second sum runs 
over all the extended states. Substitute Eq. (fT5|) into Eq. (fT^ . we have 

^i^^^(l>A, + ^ih^^^B^ = ^aA,EA,(l)A^ + ^aA, ^ f/(r - 7^p, Up)0Ai (16) 

Ai Bi Ai Ai p^Da^ 

+ ^0A,Mai + ^^-J^(yjaAr)(yj^Aj +^^-^aA,V^j(l)A^ 

Ai Ai j ^ Ai j ^ 

Bi Bi j ^ Bi j 

+ E E -|:(v.6B.)(v,fo.) + ^ -^"b. v?fB, 

Bi j ^ Bi j ^ 

In extended states {C,Bi}, we neglected the dependence of on the vibrational displacements, 
the last two terms disappear. 

Differences between the localized states and the extended states are reflected in Eq. flTBl) . 
For a localized state, we have separated 

J2uir-nn,n-)= X f^(^-^n,u")+ X t^(r-7^p,uP) (17) 

n neDA-^ P^-D^i 



the second term leads to the transfer integral which causes transition among states. The 
wave function of localized state Ai is confined in Dai ■ For a nucleus outside Da^ , its effect 
on Ai dies away with the distance between the nucleus and Da^ ■ Eq. fll7p is a generalization 
of Holstein's treatment^ to a localized state which occupies several sites. While for extended 
states we have resolved 

A/- 3A/- „ 

n n=l j=l 

the second term is the electron-phonon interaction. An extended state spreads over whole 
sample, it feels the vibrations of all the atoms. Eq. (fT8|l is similar to the usual treatment of 
e-ph interaction in a crystal. 

The first sum in Eq. fll7p only includes the atoms in the region where the wave function 
of localized state Ai is nonzero. If a localized state is close to mobility edge, it extends to a 
very large spatial region. For such a state, Eq. f|T7|) is not very different from Eq. ffTSl) . 

By applying f d^rcj)*^^ to both sides of Eq. (fT6|) . one obtains: 

d d 
^Ta^aA^^qI ~ ~ h,,}aAx ^^^mbA^^q: - - K}hBx = (19) 



Aa^i'^Ai 



Ai j ^ Ai j ^ 

Application of / d^r^*^^ to both sides of Eq. (fT6l) yields: 

d d 
E Yb^aA^^-q^ - Eax - K}aAx + {^^-qI ~ - K}bB2 (20) 

Ai 

dhCBAJ{T - 7^p, UP)0^, + ^ ^ X, / d^CB^Q^iB, 

+ EE-^(^^«^^) / ^'^^B.(V,0Aj + $^$^-^a^, [ d^CBX^^A. 

Ax j ^ Ai j ^ 

The coupled equations Eq. flT^ and Eq. fl2U]) describe the time evolution of the states un- 
der the influence of the vibrations: the transition between two localized states (LL), the 
transition between two extended states (EE), the transition from an extended state to a 
localized state (EL) and the transition from a localized state to an extended state (LE). 
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In other words, for a given initial state, i.e. {a^^(xi, ■ ■ ■ , x^j^-^ t = 0), Ai = 1, 2, • ■ ■ , A^^} 
and {hB^{xi, ■ ■ ■ ,X3_\f;t = 0),Bi = 1,2, ■■■ ,Ne}, at any subsequent moment t > the 
state of the system is completely determined (i.e {a^^ (xi, ■ ■ ■ , x^j^j-; t),Ai = 1, 2, ■ ■ ■ , N^} 
and {bB^ixi, ■ ■ ■ ,X3j\j-;t),Bi = 1,2, ■■■ ,Ne}) by Eqs. (HM and fl2Ul) . The number of vi- 
brational degrees of freedom is macroscopic, so that {aAi(^ = 0),Ai = 1,2,--- ,Nl} and 
{bsiit = 0, -Bi = 1, 2, ■ ■ ■ , Ne} cannot be assigned precisely. A description based on density 
matrix is more appropriate. 

Let us estimate the order of magnitude of the last two terms in RHS of Eq. (fT9!) or Eq. (!20!) . 
From the 1*^* order perturbation correction about the nuclear displacements to the localized 
wave function, we have 



^^'^^^ = ^ E,^ - Ec '^^ ~ ^ ^^^^ 



c 

where Ve is a typical velocity of electron, C is eigenvalue A or 5 of ha- To reach last step, 
we notice the derivative of f/(r — TZj) respect to suitable coordinate component of electron 
equals — Thus ~ PeU/h, Pe is momentum of the electron, and Ve = A typical 
term in the 2"^^^ sum from last in RHS of Eq. ( fT9l) becomes 

^V,a^,-V,0^, ~^f/ (22) 

where f„ ~ Mj^hV j is typical velocity of a nucleus. A typical term in the last sum in RHS of 
Eq. (fT9l) can be estimated from 2nd order correction to (pAi about the nuclear displacements 

V -0^, U (23) 

In the last ~, we applied the virial theorem mvl ^ U, m is the mass of electron, U is the 
interaction potential energy between electron and some nucleus. Clearly 

where a; is a typical vibrational amplitude of atoms, d is the distance between two nearest 
neighbor atoms. Notice Vn/ve ~ 10~^, m/M ~ 10~^ and x/d ~ 10~^ — 10~^, combining Eqs. 
(1221), (El and (El, we find 



$^vr(V,«Aj(V,0Aj or ^ V|0Ai « U{r-n^,M^)<PA, or ^a;,- 



(25) 



The 3rd term and the 4th term in the RHS of both Eg. (1191) and Eq. (l20p can be ignored. 
Eq. f[T9l) is reduced to 

d d 

^TA^Ariih— - Ea, - K}aA, + X]^^2Bi{^^^ ~ ^^i ~ ^J^^i = (26) 

In the right hand side (RHS), each term in the first sum is a transition between two locahzed 
states mediated by a transfer integral. Each term in the second sum is a transition from 
an extended state to a locahzed state caused by the electron-phonon interaction. Similarly 
Eq. fl20l) is reduced to 

d d 
^ Yb^aA^^-q^ - Ea, - K}aA^ + {^^g^ - - K}hB2 = (27) 

rfVe^^,f/(r-7^p,uP)0A, + 5^65,5^x, / dhCB2Q^iB, 

»1 P^UA^ Bi j 1 

In the RHS, each term in the first sum is a transition from a localized state to an extended 
state induced by a transfer integral. Each term in the second sum is a transition between 
two extended states caused by electron-phonon interaction. 

Eqs. (l26|) and (!27|) are completely general. To simplify them we need two connected 
technical assumptions (i) Yb^a^ 

~ << 1 and (ii) SA2A1 1- Obviously they are 
not true in general. For localized states which are close to mobility edge, they spread in 
many distorted spatial regions. The overlap integral YbjAi between them and an extended 
state is not small (when e-ph interaction appears, they are no longer eigenstates of ha). 
Assumption (i) means that we do not consider the localized states very close to mobility 
edge and consider only the most localized states. The contributions to conductivity from 
near-edge localized states are ignored. For the most localized states (the most low-lying ones 
in conduction band), which spread over only several bond lengths at most, condition (i) is 
satisfied. Condition (ii) is satisfied for two localized states which do not overlap. It means we 
exclude the indirect contribution to conductivity from the transitions between two localized 
states with overlapping spatial regions. For two eigenstates A2 and Ai of static hamiltonian 
ha, the overlap integral between two states SA2A1 = 0. In addition, if the distance between 
two localized states is larger than one bond length, the overlap integral can be neglected 
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even one takes into account electron-phonon interaction (e-ph). When e-ph interaction 
is taken into account, A2 and Ai are no longer eigenstates of ha, SA2A1 is not negligible 
when the spatial regions of two localized states overlap. For a localized state which spreads 
in several distorted regions^^'^^, its wave function is not single exponential decay function 
which is only suitable for the most localized statesi^. The overlap integrals YB2A1 and S^j^i 
involving such a localized state are oscillatory. To make a semi-classical estimation, one 
needs detailed information of wavefunction which could be obtained through a WKB-like 
exponential transform. We wish to remove two assumptions in a later communication^^. 

The transition between two localized states is significant only when the distance R12 
between the two is not very large. For two localized states A2 and Ai, SA2A1 << 1 if -Dai 
and do not overlap. The terms multiplied by SazAi can be neglected for localized states 
which their spatial regions do not overlap. What is more, the transfer integral is important 
only when the atoms p ^ Dai fall into Da2 or Ai = A2, 



Ax P^C'Ai M Ax 

where 

Ja2M = j d'r^X E U{r-np,nP)^A„ W^a. = j ^ ^7(r-7^p,uP) (29) 

Here, the Wa2 term only affects the self energy of a localized state through OAa- Eqs.( 12K]l and 
f l2U]) are a generalization of Eqs.(I-14) to (1-16) of ref.J. Comparing with Eaj and with 
Wa2 can be neglected. JA2A1 causes the transition from Ai to A2, it comes from the attraction 
on the electron by the atoms in Da2- For those most localized states, the wave functions 
take form of 0Ai ~ e-l^-^-^il/^i. J^aAi is estimated to he^ - ^^^J/J (l + ^)e-^'^/^, where 
average localization length ^ is defined by 2^^^ = + ,^2^^. R12 is the average distance 
between two localized states. Eg is static dielectric function, Z* is the effective nuclear charge 
of atom, Na2 is the number of atoms inside region Da2- 

Similarly JaiAz causes the transition from A2 to Ai, it comes from the attraction on 
the electron by the atoms in -Dai- For two localized states Ai and A2 in different regions, 
no simple relation exists between JA2A1 and JaiA2- This is in contrast with the situation 
JA2A1 = {JA1A2)* of small polarons in crystal where translational invariance exists^. Later 
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we neglect the dependence of JA2A1 on the displacements {x} of the atoms and consider 
JyijAi as a function of the distance R12 between two localized states, localization length ^1 
of state Ai and localization length ^2 of state A2. 
Eqs. fl2B]) and fl7r|) are then reduced to 
d d 



and 



9 d 

^Ys^AiiilT'-g^- EA^-K}aA^ + {ih—- EB.,-K}bB2 = ^ J'uaAiaAi +^ -ft'saBi^Si (31) 
Ai Ai Si 

where 

is a linear function of atomic displacements Xj. It causes type (3) transition from an extended 
state to a localized state. If we approximate extended state as plane wave ^Bi ~ e*'^^!^, 



^AaSi ~ inelLeA^ ^^~^^B^iA2) S where m ~ ^ or is typical amplitude of vibration 

at high or low temperature. So that K'j^^^^u/ Ja^Ai ~ e^^'^^^u/^. The distance between two 
nearest localized states is about several A in a-Si, K'^^^Bi"^ is several times smaller than J^aAi- 
From Eqs. fl30|) and fl3T]) . type (2) transition from a localized state in region Da^ to an 
extended state is caused by the transfer integral 

J'b2A,= Yl I d?rCB^U{r-n^,u^)4>A, (33) 

Jioi by the electron-phonon interaction K'j^^^_^ (although we will see the process does involve 
several phonons in high temperature regime in Sec. IV, it consists with the intuitive picture of 
'phonon induced delocalization'^^). J' B2A1 does not involve atomic displacements explicitly. 
Later we neglect the dependence of J'bjAi on the displacements of atoms and only view 
J' B2A1 as function of only. If we approximate extended state as plane wave ~ e~*^^2^, 
J' B2A1 ~ A_-Keoe\A ~^ ^^B2^Ai)~^- J' B2A-^ IS in the same order of magnitude as Ja2Ai- 
According to Eqs.( l30i) and ( l3Ti) . J' B2A1 does not create transitions from an extended state to a 
localized state. The asymmetries in Eq. (l32il and Eq.( |33l) come from the different separations 
Eq. (ITTI) and Eq. (fT8|) of the single particle potential energy for localized states and extended 
states. One should not confuse with the usual symmetry between transition probabilities 
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for forward process and backward process computed by the first order perturbation theory, 
where two processes are coupled by the same interaction. 

Type (4) transition between two extended states Bi and B2 is caused by electron-phonon 
interaction: 

d^rCB^gj^^Bi, Kb^Bi = {KBiB2y (34) 

It is almost the same as the usual scattering between two Bloch states in a crystal by the 
electron-phonon interaction. If we approximate extended states C,Bi and C,b2 by plane waves 

with wave vector fci and ko, KboB^ ~ 1 t-^tw— r^, where n ~ — l^i is the Thomas- 

Fermi screening wave vector. In lightly doped or intrinsic semiconductor, /t is hundreds 
even thousands times smaller than 1/a, a is bond length. Since for most localized state, 
localization length is several times a, KB2B1 ~ {'^O'^-^mbi^ much weaker than three 
other coupling constants. Contrast with type (1), type (2) and type (3) transitions, the 
transition probability of the transition Bi —>■ B2 equals to B2 Bi. two processes are 
coupled by the same interaction as illustrated in Eq. fl34|) . Eq. fl30|) and Eq. fl3T]) correspond 
to the generalized master equation for the reduced density matrix of an electron in a phonon 
batb22. 

In the conduction band of an amorphous solid, the energies of extended states are higher 
than the energies of the localized states close to the bottom of the band. If a localized state 
is not very close to the mobility edge, its localization length is small. The overlap between 
it and an extended state ^62^1 

~ may be neglected. Eq. (!30|) and Eq. (1311) read as 
Ea2 - K}aA2 = ^ JA2AiaAi + 

Y^K.bM (35) 

Ai Bi 

Eb2 - K}bB2 = ^ J' B2A^CiA^ + ^ KB2BibBi (36) 
Ai Bi 

If a localized state is close to the mobility edge, the overlaps between it and the extended 
states are not necessarily small. Eqs. (135|) and (136!) could be used generically with the 
restrictions (i) and (ii). 

In Eq.f l35l) and Eq. fl36l) . the e-ph interaction for localized state and extended state has 
been treated differently. The reason is as following. Let us consider the potential energy of 



and 



(4 
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atoms 



2^^^i'=^i^^ 9p^P^ dp— (37) 

jk p P 

the first term is the mutual interactions among atoms. The second term, the e-ph interaction, 
acts hke an external field. The first member of Eq. (1371) can be written as: 

1 . 1 - 

9 ^jki^j ~ ^'j)i^k — x^) — kjkX^xl (38) 



jk jk 



where 



^0 



J2ik-')'npgp,m = 1,2, ---SX (39) 



is the static displacement for the m^^ degree of freedom. The constant force — exerted 
by electron on the p^^ vibrational degree of freedom produces a static displacement for 
the m*^ degree of freedom. The deformation caused by the static external force of e-ph 
interaction is balanced by the elastic force. A similar result was obtained for continuum 
modeU^. The last term in Eq. fl38|) is the polarization energy, a combined contribution from 
elastic energy and e-ph interaction. The static displacement caused by e-ph interaction 
is meaningful only when the static displacement is comparable or larger than the thermal 
vibrational amplitude ■\J jf^ and the zero point vibrational amplitude \J~^^- For a localized 
state, one needs to make following substitution 

J dp if p G Da, 
g.p^dp= i (40) 

y li viDA, 

in corresponding formulae. Because the wave function of a localized state is concentrated, 
the e-ph coupling for localized states is treated as a static displacement, a;^ ~ 4^^*g2fc ; k 
is the typical value of spring constant of a bond. The e-ph interaction for extended states, 
which cause scattering among extended states or from extended state to localized state, is 
treated perturbation. 



1. Electron transfer 

If there are only two localized states in the "one electron -|- many nuclei" system, Eq. (!35 
is simplified as 

'hd^'\=Hj^A,H,^, = h,J'^^-^'--^''^' I (41) 



dt 
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Eq. (HT!) is more general than the usual spin-boson model of electron transfer. The reason 
that B.tot is not Hermitian is discussed in the paragraph adjoining Eq. fl29l) . Rtot in fHT|) can 
be separated in a different way: 

Htot = Hel + Hrc + Hj + Hb (42) 

where 

He. - I I (43) 



J21 El 



is a bare electronic two-state system, and 



- J] c/pa;p I ^ ° - J] dpXp f ° ° (44) 

describes the coupling of the two electronic states with the vibrational degrees of freedom 
(reaction coordinates) of the two localized states. The degrees of freedom which belong to 
Di or D2 are reaction coordinates, removing these degrees of freedom in h^, one obtains the 
bath Hamiltonian 

^^ = E'-^^? + ^E'w^ (45) 

j ^ jk 

The primes on the two summation signs indicate that the degrees of freedom of the atoms 
inside Di or D2 are excluded. 

1 X \ 

= 2 ^ ^ %fea;jXfe + ^J2Y1 ^ikXjXk (46) 

fceDi j^Di k&D2 j^D2 

is the coupling of reaction coordinates with the degrees of freedom of the thermal bath. 
Eq. (H2!) is the starting point of electron transfer theory^^. In this representation, the hopping 
of an electron between two localized states appears as that a particle with mass about atomic 
mass moves along reaction path between two wells, Hb is the environment in which the 
particle moves. Hi describes the diffusion of the particle and the friction suffered by the 
particle. 
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2. Capture and escape of a particle 



If there is one localized state and one extended state, Eqs. (l35i) and (136|) are simplified to 

,H,e = K+\ ^^^^^ I (47) 





1 = i^ce j 


(bs\ 




V j 







The reason why Hce is not Hermitian is discussed after Eq. (l33|) . The capture-escape Hamil- 
tonian can be resolved in a traditional way: 

H,, = H,j + H'^c + H'i + H's (48) 

where 

m, = ( ] (49) 

is a bare electronic two-state system: one is the localized state, another is the extended 
state. 

describes the coupling of the two electronic states with the vibrational degrees of freedom 
(reaction coordinates) of the localized state. The degrees of freedom which belong to D2 
are reaction coordinates, removing these degrees of freedom in hy, one obtains the bath 
Hamiltonian 

^^ = E'-^^? + ^E'w. (51) 

j 1 jk 

The primes on the two summation signs indicate that the degrees of freedom of the atoms 
inside D2 are excluded. 

= \ J2 Yl kjkXjXk (52) 

is the coupling of reaction coordinates with the degrees of freedom of the thermal bath. 

In this representation, the transition of an electron from a localized state to an extended 
state appears as that a particle with mass about atomic mass diffuses in a viscous liquid like 
a Brownian particle and escapes the trap of a well. Taking a semi-classical approximation, 
one obtains Kramers' problem of particle escape across a barrier—. The transition from an 
extended state to a localized state appears as that a well captures a particle moving in a 
viscous hquid. 
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D. Normal Coordinates 



As usual it is convenient to change the displacements {x} of the atoms to normal 
coordinates^^i^i {9}, 

Xk = J2 ^kaQa, {A^kA)po. = do^pM^ul, « = 1, 2, • • • , SAT (53) 

a 

where A^,^ is the minor of the determinant \kik — uj'^MiSikl = 0, A-^ is the transpose matrix 
of the matrix (Afc^). The two coupling constants in Eg. (1521) and Eg. (1311) which involve e-ph 
interaction are expressed as 

K'a.b. = E , KZb, = E / d'r^M (54) 

rv T 



and 



Kb,b, = E ©"^^.B, , = E / (55) 



where K'^^^^ and i^'^^^^have the dimension of force. Eg. (l35l) and Eg. (l36|) become: 
d 



d 

ih—-hA2)aA2i---Qa--- ]t) = ^ J^^^^OAi +E^^2i?i^Si (56) 



and 



where 



d 

i^^-ET - ^32)^32 (■ ■ • ©a ■ ■ ■ ; ^) = E J'B2A^a'A^ + E ^i?2Bi^Bi (57) 



= Ea^ + K, hB2 = Eb2 + (58) 

hAi describes the polarization on the amorphous network caused by an electron in localized 
state (pA^ through e-ph coupling. (" ■ ■ ©a ■ ■ ■ ; ^) is the probability amplitude at moment t 
that the electron is in localized state A2 while the vibrational state of the atoms is given by 
normal coordinates {O^, a = 1, 2, ■ ■ ■ , SA/"}. &b2(' ! ^) is the probability amplitude 

at moment t that the electron is in extended state B2 while the vibrational state of the 
atoms is given by normal coordinates {Oo, a = 1, 2, ■ ■ ■ , SA/"}. 

Due to the coupling of localized state Ai with the vibrations of atoms, by a similar 
reasoning as Eg. (1391) . the origin of each normal coordinate is shiftecl^^'^ 

©a^©a-©fS Qt' = [MctolY^ E dpA,\A,o.. (59) 
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where ~ ^ ^^^^^ fj^ ^2 \ is the static displacement in normal coordinate of the 
mode caused by the coupling with localized state Ai, where Na, is the number of atoms in 
region Da,- Eq. (!59|) leads to a modification to the phonon wave function (cf. Eq. (!6T|) ) and 
a change in total energy (cf. Eq. fl62|) ). Using {k~^)jk = {k~'^)kj and the inverse relations of 
Eq. (l53l) . one finds that the shift G^^ of origin of the a*^ normal coordinate is related to the 
static displacements by 

Qt^ = J2i^-'Uxl xIeDa, (60) 

A; 

The eigenfunctions of Ha, are 

SAT 
a=l 

where Hi^{z) is the A^*^ Hermite polynomial, 9a = { '^'^°j^°' Y^'^Qa is the dimensionless normal 
coordinate and 9^"^ = (^Mai£<i)i/2Q^i_ fjj^g corresponding eigenvalues are 

4^^ = < + E(^- + ^)^^"+^l' ^l = -^E^-^'(®:^^)' (62) 

a a 

Ea-, ~ In a small polaron, the electron is mainly localized on one site. The lo- 

calized electron deforms the crystalline lattice. In an amorphous semiconductor, the electron 
in state Ai polarizes the network, the energy of state Ai is shifted downward by Eq. (!62|) 
is an extension of Holstein's small polaron theory. Observe that S\_^ does not depend on the 
number of phonons in each mode and is a static property of the amorphous solid. £ai 
derived in usual Born-Oppenheimer approach and is called "dynamic potential energy"—. 
With the help of Eq.f l39l) and the inverse relation of Eq. fl53l) . one can show that the two 
expressions for the polarization energy are the same: 

2 ~ 2 dpdq{k ^)qp = \S\^\ (63) 

jk pq 

The continuum form of Eq.f l63l) is Eq.(6c) of refJlOl. The eigenvalues and eigenvectors of /i^i 
are 

SAT 

We change probability amplitude a^^(- ■ ■ 0a ■ ■ ■ ; t) from representation of normal coor- 
dinates to representation of occupation number, i.e. expand with the eigenfunctions of Ha, 

aA, = ;t)4f ^e-*<"'/^ (65) 



■N' 
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(7^1 (. . . A*"^ ■ ■ ■ ',t) is the probability amplitude at moment t that the electron is in local- 
ized state Ai while the vibrational state of the nuclei is characterized by occupation num- 
ber {N'^, a = 1,2, ■ ■ ■ , 3A/"} in each mode. Similarly we expand the probability amplitude 
bsii' ■ ■ ©a ■ ■ ■ ',t) with eigenfunctions of /ibi 

= E ^""^ (■ ■ ■ ■ ■ ■ ' ^e-<^*/^ (66) 

...TV' --- 

pBi . . ]\f'^ ■ ■ ■ ; t) is the probability amplitude at moment t that the electron is in extended 
state Bi while the vibrational state of the nuclei is characterized by occupation number 
{N'^, a = 1,2, ■ ■ ■ , 3Af} in each mode. 

Substitute Eq.([65]) and Eq.([66]) into Eq.([56D and applying / JJrf^^^Jj^"^ to both sides 

a 

we obtain 

+ E (^2 ■ ■ ■ A^a ■ ■ ■ \Ve7'''\Bi ■■■K--- (■ ■ ■ AT^ ■ ■ ■ ; t)e'*(4^"'-<"*)/^ 

where 

{A,...N^... \VrMi ■■■K---) = Jma, I n^^a^if ^ (68) 

describes the transition from localized state Ai with phonon distribution {■ ■ ■ N'^ - ■ ■} to 
localized state A2 with phonon distribution {■ ■ ■ ■ ■ ■} caused by transfer integral J^j^i 
defined in Eq.(l29]). 

{A,---N^---\Vlr,^^\B,...K...)= j n«i^^(E0"^%Bj4f^ (69) 

a a 

is the transition from an extended state to a localized state induced by electron-phonon 
interaction. 

Similarly from Eq. fl57p we have 

^^^= E (52■■■A^„...|F-|A•••^^^■■)Cf^,/<"^"<^^^^^ (70) 
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where 

(52 • ■ ■ iV„ . . . \Vt^\A^ ■■■K---) = J'b,a, I n^e^^f ^^if ^ (71) 

describes the transition from locahzed state | ■ ■ ■ A^^ ■ ■ ■ ) to extended state I-B2 ■ ■ ■ ■ ■ ■) 
caused by transfer integral J'b2Ai, the dependence on {xj} in J' is neglected. 

{B,...N^... \V^7'\Bi ■■■K---)= j Wde.'E.f^XY^ QaKlBM''^^ (72) 

a a 

is the matrix element of the transition between two extended states caused by electron- 
phonon interaction. It is similar to the usual expression in a metal. Eg. (1671) and Eq.( 170l) are 
the evolution equations in second-quantized form. The phonon state on the left hand side 
(LHS) can be different from that in the right hand side. In general, the occupation number 
in each mode changes when the electron changes its state. 

III. TRANSITION BETWEEN TWO LOCALIZED STATES 

If we only consider localized states, Eq. (l67|) is simplified to 

Ai-N^- 

Eq.dSHD is reduced to 

d 

{ih— - hAoJttAsi- ■ -Qa- ■ ■) JAsA^aA^ (74) 

Ai 

Eq. (1731) corresponds to the small polaron problem^'^'^ and Eq. flTl]) corresponds to the electron 
transfer among ions in polar solvent^i^i^. The transfer integral Ja^a^ decays exponentially 
with the distance between the states Ai and ^3. According to Eq. fl59|) . the shift of origin 
B^^ for the a^^ mode is order of A/" is the number of atoms in sample, A^^3 is the number 
of atoms in -D^a- If ^3 and Ai are not close to the mobility edge, (^^^ ~^aO is infinitesimal. 
To compute the integral in Eq. (l^ . we expand (^jqt^iOa — 6^^) around {9a — 0^^) to second 
order in {9^^ — 9^^). To the second order of the small parameter (^^^ — 9^^), the integral 
for the a^^ mode in Eq. fl68l) can be effected. 

A. Ja^a^ perturbation 

If the transfer integral Ja^a^ small, the semi-classical Fokker-Planck equation is 

an useful approximation^'^ for solving Eq. (l7^ . However in amorphous solids, the transfer 
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integral Eq. (l29|) between two localized states can be small. Perturbation theory can be used 
to compute the transition probability^ W from state ^^^/ y to state "^"^^ y If in the initial 
instant (t = 0) the "one electron + many nuclei" system is at state |Ai ■ ■ ■ A*"^ ■ ■ ■ ), then only 
(7^1 . . . . . ; t = 0) = 1 and other coefficients are zero. At later moment t, probability 
amplitude C^^i^A^, ■ ■ ■ ■ ■ ■ ,t) is simply determined by Eg. (175]) . The transition probability 
per unit time is 

W{A,...K...^As---N^---) = ^\{As---N^---\Vll\A,...K...)\' (75) 

1 1 dt' jl dt" exp{^(£f^^^ - 8f^,J} 

Because there exist low frequency acoustic modes in any solid, limit t —>■ (yo must be taken 
at a later staged. For two localized states, transfer integral (!29|) does not satisfy Ja2Ai = 
(-^^1^2)*- Oils must carefully distinguish the transition \Ai ■ ■ ■ N'^ ■ ■ ■) \As ■ ■ ■ ■ ■ ■) 
which is caused by the attraction on the electron by the atoms in D^^ and its inverse 
transition IA3 ■ ■ ■ Na ■ ■ ■) \Ai - ■ ■ ■ ■ ■ ) which is caused by the attraction on the electron 
by the atoms in -Dai- In the usual situation, the perturbation which leads to the forward 
transition is the same as the perturbation which leads to the backward transition. The 
two transition probabilities equal. If the spatial regions of two localized states have many 
overlaps, the difference in the transition probabilities between the forward and the backward 
direction becomes small. 

Substitute Eq. (l62l) and Eq.(|68ll into Eq. (!75l) . notice that the product of two conflicting 
Kronecker delta symbols is zero, the transition probability becomes 

W{A, ...N'^...^As---N^---) = ^ Ijt' exp{'j[{El + 8^^) - « + 8^^)]] (76) 

a 

where terms order of {9^'^ — 9^^y and higher are neglected. 

Next we sum over all possible final phonon states {■ ■ ■ A^^^ ■ ■ ■ } and take a thermal average 
over initial phonon states {■ ■ ■ A"^ ■ ■ ■ }, the transition probability from Ai to ^3 is 

WiA,^As) = Z-' Yl W...n;^...{Ai As)exp[-pJ2iK + ^)nw^] 
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exp{--J2iC^ - ^a'?[coth^{l - cost'u^) - tsint'uJa]} (77) 



where Z = XI-.-at' ■ ■ 6^P["~/^ + is the partition function of phonon. Here 

we imphcitly assumed that the vibrational state of atoms stays in equihbrium, and is not 
affected by the motion of other electrons. 

Because — 0^^)"^ ~ A/" is large and the integrand is analytic about t', the integral 

in Eq. (!77|) may be estimated by the method of steepest descent. Let the time derivative of 
the exponent of last exponential of Eg. (1771) equal zero, one has coth^^ = zcott't^Q,. The 
saddle point is t'^ = To effect the time integral in Eq. lFTTl) . change variable^ from t' to r : 
t' = ^ + T. The path of integral is also changed dt' =^ J_i_m + dr + ^ dr, 
Eq. fl771) becomes: 

W{A, As) = ^Mt) + g^iit) + 92n{t)] exp{^[{El + 8^') - {E^ + 8^^)]} (78) 



where 



giit) = £^rexp{^[(E°3 +0 - « +^^)]} 



e^P{-2 - Cricoth - csch^ cos™.)} (79) 

a 

g2i{t) and g2ii{t) are obtained by taking the upper sign and the lower sign in 

±.exp{^[« + Sf') - (El + O]} ^ ' dseM^KEl + 8^^ ' (^1 + ^^)]} 
exp{-^ J2iet^ - e)'(coth ^ - csch^ cos(t ± ^s)uJ^)} '^^^^ r'/' (80) 

a 

respectively. Since we are only interested in 'long' transition time, g2i(t) and g2ii{f) will be 
discarded hereafter-. 

If the initial phonon state is the same as the final phonon state, Eq.f l76l) is changed into 



wi^\A, ^ A,) = -i|iLexp{^|(BX, + S^') - + Ol} (81) 



rfrexp{-[(E°3 + ^,^3) _ (i^o^ +£:,^^)]}exp{--5^(e -e)'coth^} 



22 



As a counterpart of Eq. (178|l . Eq. (l8T|) is only an intermediate step in effecting the time integral 
and does not have any physical meaning. In deriving Eq. flTSl) . i.e. the steps of summing 
over all final phonon states and thermal average the initial phonon states, both {N[ ■ ■ ■ N^j^/) 
and (A^^i ■ ■ ■ N^j^f) are treated as independent variables. Subtracting Eq. flST]) from Eq. flTSl) . 
the phonon assisted transition probability from ^3 to Ai is 

Wt{A, ^ A3) = ^ exp{^[(E°3 + - iE% + S,^^)]} (82) 
exp{--$^(e-e)'coth^} 

a 

f dreMjUEl + S^'^) - (El + ^^)]}[exp{i J^i^^^' ' e)'csch^ cos™.)} - 1] 
B. High temperature limit 



0.1 




250 300 350 40( 

T (K) 



FIG. 1: (Color online) Mobility from LL transition in a-Si (T> 250K): solid line is computed from 
Eg. ([85]) with parameters: Xll = 0.2eV, AGll = 0.05ey, Jll = 0-02eV. The star symbols are 
experimental data^i^i^i^. 

Even at moderately high temperature fc^T > hZJ, the contribution from diagonal tran- 
sition, i.e. 1 in last square bracket can be neglected. The main contribution of the 'time' 
integral comes front^ the neighborhood of r = 0. Expand cosrct;. ~ 1 — and effect 

the Gaussian integral, Eq. (15^ is changed into 

Wt{A, - ^3) = %^ exp{^[(i?°3 + S^^) - « + S^^)]} (83) 
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1 PfhJa 

exp{--5^(e-e)'tanh^} 

a 

It can be viewed as a generalization of Holstein's result for small polaron hopping (cf. Eq.(77) 
of ref.=) to the transition between two localized states. The transfer integral Ja^Ai oc e~'^^'^^ , 
q is some constant. In the electric conduction, only the jumps between the nearest localized 
states are important. The dc conductivity from the localized electrons is^ 

c^dc^ = /iho>Lc(l - c)e, f,^^^ = ^WT (84) 

where ni is the number of electrons in the localized states, c = til/Nl. Eg. fl84l) was derived 
from the Kubo formula for optical-phonon-assisted hopping^. 

At 'very' high temperature ksT > 2.5^cJ, one can use tanhx ~ x and cschx ~ ^ (the 
error is less than 0.003 when x < 0.2). Eq. (!83l) is changed into 

WTiA-.A,) = u,,e-^^'/'^^, -ll = ^[y^Y^\ El:' = ^{l + ^f (85) 

h ALLksT 4 All 

where 

AGL = (^^3 + ^b^') - (El + ^^), >^ll = IY1 MaU^Uei^ - Qt^f (86) 

a 

AG^LL is the energy difference between two localized states. All is the reorganization 
energy which depends on the vibrational configurations {6^^} and {O^^} of the two localized 
states. Because 0^ does not have determined sign for different state and mode, one can 
only roughly estimate All ~ ^"^(if^;^)^- i^ about several tens meV. It is in agreement 
with observed value of a-Si. According to Eqs.l HOi) . ( l39l) and ( 160|) . if two localized states are 
in two different spatial regions, then the atoms displaced by e-ph interaction are different 
for two states. All can not be zero. If the spatial regions for two localized states have 
overlap, there is possibility that some atoms are moved in similar way in two states by e-ph 
interaction. All could be small. All is zero for two extended states, because for extended 
states the static displacements of the atoms are zero. Molecular dynamics gives phonon 
spectrum and the eigenvector Ap^^Q, of each mode, and therefore the reorganization energy 
All- Eq. (!85|) is the transition probability between two localized states in an amorphous 
solid, it has the same form as Eq.([T]). 
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Fig{T]is the temperature dependence of mobility of a-Si for T> 250K(w0.02eV), plotted 
from Eqs. flH^ and fl85|) according to a group of typical parameters: Ra^Ai = SA, Xll= 0.2eV, 

AGll = 0.05eV (about one half mobility edge) and Ja^Ai = 0.02eV. Using above parame- 
ters, Wt ~ lO^^sec"^. The computed mobility agrees with the measurementa^'^'^i^. 

C. Low temperatures 

At low temperature, the probability of diagonal transition and the probability of non- 
diagonal transition is comparable. The probability of diagonal transition is given by 

W^'\A...N^...^A,...N^...) = ^Jl,AAiEl + - {EX + ^^)] (87) 

exp{-5^(iV. + i)(e-e)n 

Q 

in the same approximation as Eq.( |76l) . The first line is the expression for the transition 
between two electronic states without the phonon environment. The last exponential factor 
of Eq. fl87|) also appears in the corresponding expression for non-diagonal transition. It can 
be explained as the decrease of transfer integral caused by electron-vibration coupling. 

At zero temperature = 0, and only zero point vibration remains. The transition 
between two localized states is then pure quantum tunneling. At finite temperature, one 
needs to average Eq.llHTl) over the equilibrium phonon distribution: 

W?{A, A,) = ^JUAiE^' + £^^) - {E^^ + £^n] (88) 

a 

The lower value of diagonal transition probability in non-zero temperature comparing with 
that of zero temperature may be understood as that thermal vibrations disturb quantum 
tunneling. Group the exponential function and JasAi -^1- dSH]) together, one may say that 
transfer integral is reduced in a multi-phonon process. This is similar to that the effective 
band width decreases with e-ph coupling which was derived from a different model^i^. 

In low temperature ksT < hZJ/lO (cschx < 0.0135 when x > 5), the exponent in the 
last exponential of Eq. (|82l) is small. The exponential can be expanded in Taylor series of its 
exponent, then the 'time' integral can be finished. Denote 

/M = ^(e-e)^csch^ (89) 
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FIG. 2: (Color online) Mobility from LL transition in a-Si (T< 250K): solid line is computed from 
Eq.(l90]) with parameters: AGll = 0.05ey, Jll = 0.02eF, w = 27r x 850cm-^ The star symbols 
are experimental data^i. 

the result is 

Wt{A, - A,) = ^^exp{^^^}exp{-l$:(e -^.^^)^coth^} (90) 

a 

a 

aa' 
+ ■■■} 

In fact coth in the second exponential factor can be replaced by 1 when fc^T < fioJ/8 
(cothx — 1 < 0.0007 when x > 4). The terms in the curly bracket are single-phonon 
transition, 2-phonon transition, etc. / quickly decreases with decreasing temperature, only 
single-phonon processes are important at low temperature. Because low frequency modes 
have larger density of states in amorphous solids, multi-phonon processes are also visible 
at low temperature. FigI2] is the mobility from LL transition at low temperature regime 
calculated from Eq.( [90i) with parameters AGll = O.OSeV, Jll = 0.02ey, = 27r x 850cm"^ 



(taken from refiSTI). No reorganization is needed at low temperature regime, but function f 
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decreases with temperature rapidly. The separation temperature T=250K in data fitting^i^^ 
is higher than the estimated upper hmit of low temperature regime tu/lO ~ 122K. 

The activation energy for the transition between two mid-gap states is just the energy 
difference between the two levels. For mid-gap states, variable range hopping is an effective 
conduction mechanism at low temerature^i^. Eq. fl90p is the conduction of tail states at low 

AG" 

temperature due to non-diagonal LL transition. The activation energy — 7^ is one half the 
energy difference between two localized states. This is contrast with the high temperature 
partner Eq. (!85|) where activation energy involves reorganization energy. Therefore variable 
range transition is also an effective conduction mechanism for tail states as non-diagonal 
transition in the low temperature regime (fc^T < fiwjVS). Thus is a rough upper limit 
temperature of variable range hopping. In a-Si, cJ/10 ~ 122K roughly agrees with the upper 
limit temperature^ of variable range hopping (lOOK). 



IV. TRANSITION FROM A LOCALIZED STATE TO AN EXTENDED STATE 

If at the initial moment (t = 0) the "one electron + many nuclei" system is at state 
1^1 ■ ■ ■ ^'a ' ' ')i then C^^ (■ ■ ■ A*"^ ■ ■ ■ ; ^ = 0) = 1 and other coefficients vanish. At later 
moment t, the probability amplitude F^'^{- ■ ■ Na - ■ ■ ]t) is simply determined by the first 
order perturbation approximation of Eg. ( 1701) : 

= {B,...N^... \Vt,\A, ..-K... )Cf;,,/^<^'-<^'^'^ (91) 

It describes the transition from localized state \Ai- ■ ■ N'^- ■ ■) to extended state 
\B2 - ■ ■ Na - ■ ■) caused by the transfer integral J'b2Ai- Expand ^N^{9a — 0^^) around 9 a 
to the second order of 9^^, the integral in matrix element fl7T]) can be effected. Then the 
transition probability from localized state Ai to extended state B2 is 

WAA, B2) = ^^e-^(^--^^.-^^.)/%xp{-l5^(e)^coth^} 

a 

£ dreMjiEs, - El - 4j}[exp{^ Y.^^^'^'^'^ch^ cosu^r} - 1] (92) 
For moderately high temperature ksT > huJ, the 1 in the last square bracket can be 
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neglected. Expand coscUq-t ^ 1 — ^-y^, the integral is a Gaussian integral 

J'2 



W^^ = :^-^e-'^(^-^-^-i-^-i)/%xp{-i5^(e)'tanh^} 



a 



For 'very' high temperature ksT > 2.5/iaJ, Eq. (!93|) is further simplified as: 

W^'' = ULEe-''""/'-^, i?f^ = ^(l + ^i^f (94) 

where 

AGlj, = Es,-{El+S'^J (95) 
is the energy difference between extended state B2 and localized state Ai. 



-LE = ^^(^-^)^/^ Xle = lj2M^u;l{e^^r (96) 



Xle is the reorganization energy for transition from localized state Ai to extended state B2. 
It is interesting to notice that activation energy E^^ for LE transition can be obtained by 
assume 6;^^ = in Xll. Transition from a localized state to an extended state corresponds 
to that a particle escapes a barrier along reaction patb^. 

From their expressions (!86l) and (196!) . we know X^e is same order of magnitude as Xll- 
AG°^ is order of mobility edge, several times larger than AG°l- J' B2A1 is several times 
larger than J' AzAi ■ The spatial displacement of the electron in a LE transition is about the 
linear size of the localized state. In general the LE transition probability is smaller than 
that of the LL transition. From Eqs.( |Mll and (l85l) . W'^^ becomes comparable with 
only when temperature is higher than k]^^{E^^ — E^^)[2\n{JLE / Jll)]~^ ■ The mobility edge 
of a-Si is about O.leV, therefore LL transition is dominant in intrinsic a-Si below 580K. 
However if higher localized states close to the mobility edge are occupied due to doping, 
there exist some extended states which satisfy AG\e ~ -^^^ these LE transitions, 

is comparable to E^^ . The LE transition probability is about 10 times larger than that 
of LL transition. For these higher localized states, using parameters given for LL transition 
in a-Si, W^^ ~ lO^^Sec"^. Comparing Eqs. llH^ and ( |96l) . from the value of hopping mobility 
0.05 cm^Sec^^V^^, we expect the mobility from LE transition is several tenths cm^Sec~^V~^ 
for a-Si, the same order of magnitude as observed 'drift mobility— 1^. 
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The probability of diagonal transition is given by 

w!f\B2---N^ . ■ ■ ■ iV„ . . . ) = '^J'\,aAEb, - {E% + 4J] (97) 

exp{-5^(iV, + i)(e)n 

a 

The last exponential factor can be explained as the decrease of transition integral caused 
by electron- vibration coupling. At zero temperature Na = 0, the transition from a localized 
state to an extended state is then pure quantum tunneling. At finite temperature, one needs 
to average Eq. (in7|) over the equilibrium phonon distribution: 

W!,'\A, B,) = yJ%.A.exp{-l5^(e)^coth^}5[EB. - {E^ + £^^)] (98) 

a 

In low temperature /c^T < fuJ /lO, the exponent in the last exponential of Eq. fl92|) is 
small. The exponential can be expanded in Taylor series of its exponent, then the 'time' 
integral can be finished. Denote 

JleM = ^(e)^csch^ (99) 

the result is 

WriA^ B,) = ^-^^^^ exp{^^^} exp{-i ^^(e)^ coth ^} (100) 

a 

{J2 fLEMliSiAGl^ + nw^) + 5(AG°E - hwo.)] 

a 

aa' 
+ ■■■} 

Formally type (2) transition is quite similar to type (1) transition. To obtain the former, 
one needs to make substitutions: J AzA^ J' B2A1 ©a^ — ©a^ ©a^ and [(-E^^ + S^^) — 
(E^3 + gAzy^ ^ \^Eb2 - {E\^ + 8^^)] in corresponding expressions of the former. However 
the physical meaning of the two are completely different: they are reflected in the transfer 
integrals Eqs.(l29]) and 
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Consider a localized state and an extended state, both of them are close to the mobility 
edge. Then AG^s is small. For such a localized state, the displacements of the atoms 
induced by the polarization of the localized state (cf. 9;^^ in Eg. (13^ ) are small. Its atomic 
configuration is similar to that of an extended state. Therefore the reorganization energy 
\le is small. The transition probability (1511) can be large. Similar picture has been suggested 
long time ago in name of phonon-induced delocalization^^. The inelastic process makes the 
concept of localization meaningless^i^ for the states close to the mobility edge. 

V. TRANSITION FROM AN EXTENDED STATE TO A LOCALIZED STATE 

If at the initial moment the electron+nuclei system is at state \Bi ■ ■ ■ N'^ ■ ■ ■) , then 
F^^{- ■ ■ ■ ■ ■ ]t = 0) = 1 and other coefficients are zero. At later moment t, probability 
amplitude C'^^(- ■ ■ Na ■ ■ ■ ;t) is determined by first order perturbation theory from Eq.( l67Il : 

= (^2 ■ ■ ■ iV„ . . . \Vf^^'\B, ■ ■ ■ K . . . ■ ■ K ■ ■ ■ ; t)e^*(4:'^'-<"')/^ (IQI) 

Expand ^N^i^^a — 0^'^) around 9a to second order of 6*^^, effect the multiple integral 
over normal modes, the matrix element fl69l) can be calculated. The diagonal transition 
probability (in which the phonon state invariant before and after the change in electronic 
state) induced by electron-phonon interaction is 

a a' 

(102) 

The probability of the transition from extended state \Bi) to localized state IA2) is 

Wt{B, ^ A,) = ^exp{-|(E°^ -i?^J}exp{-i5^(e)'coth^}(/i + J2) 

(103) 

where 

h = f drexp{j{El + - E^J} exp{i J](e)'csch^ cos^„r} (104) 

a' 

a' 
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and 

(105) 

2 

At moderately high temperature fc^T > /SJ, by expanding cos a; « 1 — the integrals 
1 1 and I2 can be carried out. One has 

W^\B,^A,) = ^exp{-^(i?°, +4, -EBj}exp{-l J](e)'tanh^} (106) 

a 

a' a' 

[Al^^J^. ) csch 2 ) (E„/^^(^«^^)^csch^)2J 

a' 

It is interesting to notice that the first two lines of Eq. fll06p have the same structure as 
Eq. fl83l) and Eq. fl93|) . The non-radiative lifetime of an electron in an extended state is 
defined as tufe = 1/ Za, ^t""- 

For very high temperature ksT > 2.5^ZI;, Eq. (11061) becomes 

Wel = UELC-^^"/'''^ (107) 

with 

j,^, - lf^^)i/2ri VfiT'"' e-"^)^l Vf ^^^"'^^^'^ coch ^^"' I ^^^"^iQ"'^\ cch^^^"'l 

a' a' 

(108) 

[2fV^-^f6'^n^cGch^^"r^ 1 (-^^2 + 4.--^^,)^ 

[2(l^a;J^^„ ) csch ^ ) (E>^^(^«^^)^csch^)2J ^ 

4M„/ 2 8 2 4 2 



and 

(109) 
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\el is the same order of magnitude as '^G'el < is smaller than -E'f^. 

^MBi^t^ is the same order magnitude as Ja^Ai, EL transition probability is larger than 
that of LL transition. With the parameters for LL transition in a-Si, W^^ ~ 10^'^ — lO^^sec"^. 
In a-Si:H, it leads to a mobility about several tenths of cm^sec~^V~^, a value between the 
observed 'hopping mobility' and 'drift mobility—*^. However the population in extended 
states is meaningful only when temperature is high enough or doping is heavy enough. 
For intrinsic semiconductor, EL transition is not important for conduction at not too high 
temperature. 

AG^2. < Las a deep consequence. From the expression ( l62l) for and the order of 
magnitude of mobility edge^, the energy difference AG^^ < is order of several tenths 
eV. For the extended states with \AG^^]^\/ Xel < 1, we are in normal regime: the higher 
energy of an extended state (i.e. more negative but still \AG%^\/ Xel < 1), the 

smaller the activation energy E^^. The higher extended state has shorter lifetime (0 < 
\^Gel\/^bT < 5 cf. FiglHl one takes XEL/ksT = 5), the electron has less time moves 
in extended state, it will contribute less to electric conductivity. For extended states with 
energies far above the mobility edge (such that \ AG%i^\/ Xel > 1), we are in Marcus inverted 
regime (lAG^^I/Z^^T > 5 cf. Figj3]): the higher the energy of an extended state, the larger 
the activation energy. The higher extended states have long lifetime (as usual experience: 
the higher energy of the initial state, the faster of the decay), and they will contribute 

AG" 

conductivity more. In the middle of the two regimes, ~ — 1. For these extended states, 

no activation energy is needed for the transition to localized state. They will quickly decay 
to the localized states. Reflected in photoluminescence, the luminescence insisting time as a 
function of luminescence frequency decreases with frequency. In a crystal, phonon assisted 
non-radiative transition is slowed down by the energy-momentum conservation law, the 
lifetime of an excited state is long. This phenomenon was found some times ago: in c-Si/Si02 
quantum well structure, the photoluminescence lifetime is about 1ms, the photoluminescence 
life time is not sensitive to the monitored wavelength^. Comparing with crystalline quantum 
well, the photoluminescence lifetime of a-Si/Si02 structure becomes shorter with a decrease 
in monitored wavelength^^ 13ns at 550nm and 143ns at 750nm. The trend is consisten with 
the left half part of Figl3l 

The preceding discussion is also applicable to the transitions from the extended states in 
conduction band to the unoccupied localized states in valence band. 
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FIG. 3: (Color online) Schematic illustration of lifetimes of states as function of energy: energy 
unit is k^T. A = 5, E/, is taken as energy zero, the vertical axis is the relative lifetimes of states 
relative to the shortest one. Marcus inverted region: E > 8, lifetime increases with energy. The 
short lifetime belt exists in 3 < E' < 8. The vertical axis is scaled by vel 

At low temperature ksT < one can expand the exponentials in Eqs. (11040 and 

( 11051) into power series. Then the integrals can be carried out term by term. To 2-phonon 
processes, the transition probabihty from extended state \Bi) to locahzed state IA2) is 

Wt{B, A,) = ^exp{-^« - i?Bj}exp[-l J](e)'coth^] (110) 

a 



2Ma'0i!a' 

X [6{E% + 4^ - Eb, + hwa') + 6{E% + 4^ - Eb, - huj^')] 



a' a" 

4 2M«/tu„/ 4 2^2 ^ ' 



X [(5(E^^ + - Eb, + huja' + huj^n) + 6{E%^ + - Eb, + hujo.' - huja.") 

+5(E°^ + - Eb, - hWa' + f^a") + 5(^1 + ^2 - - - ^"")] 

64 ^ 2 
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X [5{E\^ + - Eb, + hWa" + hWa»>) + 5{E\,^ + - Eb, + fkUa" - hiOa'") 

+ - Eb, - huJa" + huJa'") + 5(^1 + " ^B, " ^a" " nWa'")] 
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a' 

where 



IelM = ^(e)^CSch^ (111) 



In conduction band, the energy of any locahzed state is lower than that of any extended state. 
Thus zero-phonon process cannot conserve energy and therefore is not possible. Because 
AG^^ < 0, factor exp[| AGg^|/(2/cBT)] increases with decreasing temperature. On the 
other hand, other factors in Eg. (11 101) decrease with lowering temperature. There exists an 
optimal temperature T^,, at which transition probability is maximum. In non-equilibrium 
phenomenon like luminescence, T^, may have some traces. 

VI. TRANSITION BETWEEN EXTENDED STATES CAUSED BY ELECTRON- 
PHONON INTERACTION 

If in the initial moment the "one electron+ many nuclei" system is at state \Bi - ■ ■ N'^ ■■■)■, 
then only F^^{- ■ ■ ■ ■ ■ ; t = 0) = 1 and other coefficients are zero. At later moment t, 
probability amplitude F^^{- ■ ■ ■ ■ ■ ',t) is determined by the first order perturbation theory 
based on Eq. flTOj) : 



= {B,...N^... Wr/'IB, . . . K . . . )Ff^, ^e^*<"'-<^ (112) 

The matrix element (172]) an be calculated from the recursion relation for z^]\f{z). The 
transition probability from extended state \Bi) to extended state I-B2) is 



W{B, B,) = '^-^J2(<Bj'j^^[^S{EB,-EB-nnj^^^^ 

a' ^ 

(113) 

where iV^ = (e^'^-' - 1)"^ is the average phonon number in the a'th mode. In a crystal, 
Eq. (lll3p arises from inelastic scattering with phonons. Although in amorphous solids one 
cannot classify the electronic states and the phonon states with wave vectors, Eq. (lll3l) is 
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similar to its crystalline counterpart of plionon inelastic scattering. The relaxation time 
t{Bi) of state \Bi) can be defined as 

= J2 W{Bl - B2)[l - (114) 

a^E ^ ^^e/i^^, /i^^ = — (115) 

m 

where n^; is the number of electrons in extended states, /i^^ is the mobility of the electrons in 
extended states. From the observed valu^^ of /idrift = 0.61cm^V~^Sec~^ of a-Si, one deduces 
relaxation time r ~ 0.35fs. Take effective nuclear charge Z* of Si as 4, bond length b =2.35A, 
the coupling constant will he K ^ iJomy^ ~ ^ lO^^N. W ~ ~ 6.7 x 10i2sg(.-i_ 

The relaxation time r = W^^ ~ O.lSps. The electron moving in extended states scattered 
by phonons will give a mobility 3 order of magnitude larger than the observed value of 
'drift mobility'. Therefore the drift of electron in extended states is not responsible for the 
observed 'drift mobility'. 

A. Four types of transitions and conduction mechanisms 

Some features of 4 types of transitions are summarized in TABLE [H With the basic 
transition probabilities for four types transitions in Sec.III-Sec.VI, the transport properties 
could be derived from various master equations (Fokker-Planck equation, Boltzmann equa- 
tion etc.), we will not further pursue this well-known but complicated procedure and only 
briefly discuss the conduction mechanisms. 

In thermal equilibrium at temperature T, the ratio R of the number of electrons in the 
extended states to the number of electrons in localized states is of order ^-^/^bT ^ p jg ^j^g 
distance between the lower mobility edge and the bottom of conduction band (for a-Si, F is 
about^ 0.2eV. At T=300K, i? ~ 4 x 10~^). However the mobility for a carrier in extended 
states is much higher than the mobility in localized states. For the extended states outside 
the short lifetime belt, the life time tufe of a carrier in one of these states is long enough. 
Under the influence of external electric field, the carrier drifts in these extended states. Thus 
at higher temperature, transition from localized states to extended states (LE) will lead to 
an appreciable contribution to conductivity af^ = , where m is mass of electron, 

ue is the number of electrons in extended states, r ~ [iy(i?i —>■ -82)]^^ is the mean free 
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TABLE I: Some features of 4 types of transitions 



transitions 


origin 


phonons needed 


activated 


role in conduction 




L^L 


J 


multi 


yes 


direct 


1012 


L^E 


J' 


multi 


yes 


direct+indirect 


IQlO _ 1013 


E^L 


K' 


multi 


yes 


direct+indirect 


1013 - 1014 


E^E 


K 


single 


no 


reduce 


1013 



time in the extended states determined by Eg. (11131) . When tufc >> r, cr^^ is reduced to the 
usual formula Te^nE/m. 

In present perturbation treatment, the total conduction comes from four processes: (1) 
the hopping among localized states; (2) the hoping from a localized state to an extended 
state; (3) the hopping from an extended state to a localized state and (4) the drift in 
extended state scattered by phonon. The carriers are pumped to the extended states by 
type (2) transition. 

The last column gives the order of magnitude of the transition probability estimated 
for a-Si at T=300K. The value of LL transition is for the transition between two nearest 
neighbors. For an intrinsic or a lightly n-doped semiconductor at not too high temperature 
(T<580K for a-Si), only the lower part of conduction tail is occupied. AGle is large, LE 
transition probability is about two order of magnitude smaller than that of LL transition. 
For an intrinsic semiconductor at higher temperature or a doped material, AG^^ becomes 
comparable with AG\i, LE transition probability is about ten times larger than that of LL 
transition. The first three types of transitions increase the mobility of an electron, whereas 
type (4) transition decreases mobility of an electron. Although a stricter solution of Eq. fll4p 
could be obtained by considering the change in the functional form of vibrational states by 
the motion of the electrons and the change in electronic wave function in the same time. 
However for a localized state in amorphous semiconductor, the time spent from one localized 
state to another (10~i°m/10^m-sec~i ~ 10~i^ — 10~i^sec) is much shorter than the stay time 
on one localized state (10~i2sec). In addition the slow multi-phonon reorganization process 
of a configuration reduces the impulse of an electron on the network when the electron 
moves. Thus one can neglect the recoil of lattice by moving an electron. This differs from 
an ionic crystal where the deformation of the lattice closely stalks the motion of the slow 
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electron. 

B. Long time and higher order processes 

The perturbation treatments of the fundamental processes given in Sec.IV-Sec.VII are 
only suitable for short times, in which the probability amplitude of the final state is small. 
Starting from a localized state we only have L ^ E process and L L process. Starting 
from an extended state, we only have E L process and E —>■ E process. For long time 
period, higher order processes appear. For example L—>-E^L^L—>-E—>-L—>-E^ 
E L etc. Those high order processes are what occur in the amorphous solids. However, 
the perturbation picture could be used to compute transport properties. In semiconductors, 
the number of carriers are smaller than the number of avialable states. The motions of 
the individual carriers can be viewed as independent. On the other hand, for macroscopic 
sample, each type of states are highly populated if we consider the whole sample. Thus if 
we concern the collective behavior of all carriers rather than one individual carrier in a long 
time period, the picture of the four types of transitions works well statistically. 

VII. SUMMARY 

For amorphous solids, following Holstein's work on small polaron, we established the evo- 
lution equations for localized states and extended states in the presence of lattice vibrations 
(Eqs. (l30l) and ( !3TI) ). To simplify the evolution equations, for localized states, we restrict our- 
selves to the most localized ones. These localized states are attached to the most distorted 
regions, and are spatially well separated. One can neglect the overlap integral between two 
of them and the overlap integral between an extended state and one of the most localized 
states. For any process involving localized states close to mobility edge, the conclusions 
obtained from the simplified evolution equations (Eqs. fl35| ) and fl36l) ) are only qualitatively 
correct. 

For short times (t << j, jr, Jy), perturbation theory can be used to solve Eqs.( l35|l 
and (|36l) . One obtains the transition probabilities of LL, LE, EL and EE transitions. In high 
temperature {ksT > hw) and low temperature (/i;bT < /jU/lO) limits, the 'time' integral in 
transition probabilities of LL, LE and EL transitions can be carried out analytically. The 
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relative errors are less than 10~^. In amorphous semiconductors (e.g. a-Si), the transfer 
integrals J and J', and the e-ph interactions K' and K are sufficiently small such that the 
perturbation theory is applicable for a meaningful time period. Although in a sample the 
motion of a single carrier is a long time and high order process. Because conduction is a 
collective behavior of many carriers, the results of perturbation theory can be used. 

The hopping motion of electron along the direction of electric field during LL, LE and 
EL transitions directly contributes to conductivity. In external electric field, an electron 
in extended states moves along the direction of electric field. EE transition deviates the 
direction of drift which is along the direction of field and reduces the transport of charges. 
LE transition increases the number of electrons in extended states, whereas EL transition 
decreases the number of electrons in extended states. They also affect conductivity indirectly. 

At 'very' high temperature ( ksT > 2.5hnj), the transition probabilities for the transitions 
of type (1), type (2) and type (3) are reduced to an Eq.([T])-like formula. It is universal for 
the processes involving localized state. If a localized state close to the bottom of conduction 
band, AGle — AGll is about the mobility edge, W^^ is one or two order of magnitude 
smaller than W^^. If a localized state close to mobility edge, AGle ~ AGll, because J' is 
several times larger than J, W^^ could be one order of magnitude larger than W^^. The 
probability of type (3) is one or two order of magnitude larger than that of LL transition. 
The reason is AG EL < 0, E^^ is smaller than while K'Q is same order of magnitude 
as J. The probability of EE transition is about 10^ times larger than W^^ (cf. Table 
I). Type (4) transition deviates the direction of the drift motion along the direction of 
electric field and reduces conductivity. This is contrast with LL, LE and EL transition. The 
relative contribution to conductivity of four types of transition also highly depends on the 
number of carriers in extended states and the number of carriers in localized states, which 
are determined by the extent of doping and temperature. 

At low temperature ( < hZJ/lO), the non-diagonal transition is still multi-phonon 
activated process whatever it is LL, LE or EL transition. The activation energy is just half 
of the energy difference between final state and initial state (AGll/2, AGleI'^-, AGelI'^)- 
The reorganization of vibration configuration is not needed, this is similar to the hopping 
between mid-gap states^^'^?, which does not need reorganization of the static displacements^. 
Thus we expect variable range hopping also works for the LL transitions which are transitions 
between tail states. 
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In a-Si for both high and low temperature regimes, the computed temperature depen- 
dence of mobihty of LL transition agrees with the observed^i^i^i^i 'hopping mobihty'. EE 
transition offers a mobihty which is three times larger than the observed drift mobility. How- 
ever, in an intrinsic or a lightly doped system, the occupation probability of the extended 
states is small. The observed 'drift mobility' may come from LE transition. 

There exists a short lifetime belt of the extended states inside conduction band or va- 
lence band. These states favor non-radiation transition by emitting several phonons, have 
smaller contributions to dc conductivity than other extended states. The fact that photo- 
luminescence lifetime in a-Si/Si02 decreases with frequency confirms the existence of short 
lifetime belt in conduction band. The phonon-assisted non-radiation transition dissipates 
the population of the excited states in amorphous semiconductors. 

Some important points are not addressed in this work. The overlap integral between 
two less localized states which spread in several distorted regions and the overlap integral 
between one less localized state and an extended state require the wave function of a less 
localized state which cannot be described by a single exponential decay function. The 
coherence between states is subtle in Eqs.f l56l) and ( 1571) . The dependence of conduction on 
the strength of external field also requires including electric field in the zeroth order in the 
evolution equations. The drag effect^ of water on the conduction of solvated DNA is not 
included in Eq.f ll4l) . We wish to discuss them in near future. 
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